Anomalous Josephson effect in p-wave dirty junctions 
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The Josephson effect in p-wave superconductor / diffusive normal metal / p-wave superconductor 
junctions is studied theoretically. Amplitudes of Josephson currents are several orders of magni- 
tude larger than those in s-wave junctions. Current-phase (J-<p) relations in low temperatures are 
close to those in ballistic junctions such as J oc sin(<^/2) and J oc <p even in the presence of random 
impurity potentials. A cooperative effect between the midgap Andreev resonant states and the prox- 
imity effect causes such anomalous properties and is a character of the spin-triplet superconductor 
junctions. 
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The internal 7r-phase shift (sign change) of pair po- 
tentials is essential for unconventional superconductiv- 
ity and is the source of the midgap Andreev resonant 
state (MARS)****. It is now known that the MARS is 
responsible for anomalous transport properties in super- 
conducting junctions*. In normal metal/superconductor 
junctions, transport properties are affected also by the 
proximity effect which is interpreted in terms of diffusion 
of Cooper pairs into normal metals. In what follows, 
we assume that normal metals are in the diffusive trans- 
port regime due to impurity scatterings. Recent theo- 
retical studies have revealed sensitivity of the proximity 
effect to the internal phase of pair potentials*!*. In nor- 
mal metals attached to unconventional superconductors, 
Cooper pairs have a sign degree of freedom reflecting the 
7r-phase shift of pair potentials. Suppression of the prox- 
imity effect is usually expected because wave function of a 
Cooper pair originated from the positive part of pair po- 
tentials cancel that originated from the negative part*&. 
Two of us, however, discussed anomalous enhancement of 
the zero-bias tunneling conductance due to the proximity 
effect in a presence of the MARS**V 

In superconductor / nomal metal / superconductor 
junctions, another phase degree of freedom affects quan- 
tum transport. Namely, the external phase difference 
across the junctions tp drives Josephson currents. An im- 
portance of studying the Josephson effect is growing these 
days because quantum interference devices consisting of 
Josephson junctions can be basis of future technologies. 
In fact, a recent experiment has tried to apply high-T c 
superconductors to coherent devices 9 . In unconventional 
junctions, the MARS is considered to have the phase de- 
gree of freedom. When MARS's are formed at the two 
junction interfaces, the external phase may modify inter- 
ference effects between the two MARS's and Josephson 
currents. The research in this direction can shed new 
light on quantum transport in unconventional supercon- 
ductors. 

In this paper, we theoretically study Josephson cur- 
rents between two p-wave superconductors through nor- 



mal metals by solving the Bogoliubov-de Gennes equa- 
tion using the recursive Green function method****. We 
show that amplitudes of Josephson currents in the p- 
wave junctions are much larger than those in the s-wave 
junctions when transmission probabilities of junction in- 
terfaces are small. The local density of states in nor- 
mal metals has a zero-energy peak reflecting anomalous 
diffusion of the MARS's into a normal metal and that 
spatial profiles of the zero-energy peak depend strongly 
on ip. As a consequence, current-phase {J-f) relations 
remarkably deviate from the sinusoidal function in low 
temperatures and are close to those in ballistic junctions 
such as J oc sin(</j/2) and J oc ((3 i2 4 iiii i i-. The reso- 
nant tunneling through the MARS in normal metals is 
responsible for such unusual Josephson effect. The ob- 
tained results imply high potentials of spin-triplet super- 
conducting junctions as coherent devices. 

We consider three pairing symmetries on two dimen- 
sional superconductors: (i) = _Ao for s-wave, (ii) 
Ag2k x k y for d xy -wa,ve, and (iii) Aok x for p^-wave sym- 
metries. Here Ao is the maximum amplitude of pair 
potentials at the zero temperature, k x — k x /kp and 
kg = ky/kF are normalized wave numbers on the Fermi 
surface in the x and y directions, respectively. Joseph- 
son currents are parallel to the x direction and junc- 
tion interfaces are parallel to the y direction as shown 
in Fig.^a). The pair potentials in momentum space are 
illustrated in Fig.QJb). An interference of a quasiparticle 
enables formation of the MARS at a junction interface 
when Afc^fc A-fc^fc < 4 ' 16 . The pair potentials in the 
d xy - and p^-symmetries satisfy the relation for all wave 
numbers. The absence of the proximity effect in normal 
metals is described by a relation*!* Ak x ,k v = —Ak x -k y - 
The pair potential in the d^-symmetry satisfies the re- 
lation. Thus the proximity effect is expected in both the 
s- and p^-wave symmetries. In Fig.^b), we classify the 
pairing symmetries into three groups by the presence (o) 
or absence (x) of the two interference effects***. Within 
p-wave symmetries, we pay special attention to the p x - 
wave symmetry because the proximity effect and MARS 
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FIG. 1: A schematic figure of a Josephson junction on the 
tight-binding lattice is shown in (a). In (b), we illustrate the 
pair potentials in momentum space, where open circles repre- 
sent the Fermi surface. The pair potentials are classified into 
three groups by the presence or absence of the two interfer- 
ence effects. 



are present at the same time. On the other hand in the 
Py-w&ve symmetry, neither is present. 

Let us consider Josephson junctions on the two- 
dimensional tight-binding model as shown in Fig. ^a). 
A vector r = jx + my points a lattice site, where x and 
y are the unit vectors in the x and y directions, respec- 
tively. The junction consists of three regions: a normal 
metal (i.e., 1 < j ' < Ljv) and two superconductors (i.e., 
~°° < j ' < and + 1 < j < oo). In the y direc- 
tion, the number of lattice sites is W and we assume 
the periodic boundary condition. Electronic states in su- 
perconducting junctions are described by the mean-field 
Hamiltonian 
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among nearest neighbor sites. We assume that t and the 
Fermi energy /i are common in superconductors and a 
normal metal. In a normal metal, on-site potentials are 
given randomly in the range of — Vj/2 < e r < Vj/2. We 
introduce insulating barriers at j — 1 and Ln , where e r is 
given by Vb ■ Two superconductors in which e r are taken 
to be zero are identical to each other. In the p^-wave 
symmetry, a spin vector of Cooper pairs d points the z 
direction. The arguments below do not depend on direc- 
tions of d. The Hamiltonian is diagonalized by the Bo- 
goliubov transformation and the Bogoliubov-de Gcnncs 
equation is numerically solved by the recursive Green 
function methocUS*ii. Josephson currents are given by 
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with r' = r + x, where G Un is the Green function and 
u> n = (2n + 1)ttT is the Matsubara frequency with n 
and T being an integer and a temperature, respectively. 
In Eq. JSJ, Tr means the trace in the Nambu space 
and the summation over m. In this paper, the unit 
of h = ks = 1 is used with ks being the Boltzmann 
constant. Local density of states is also calculated from 
N(E,j) = — Im Tr G,E+i 7 (r, r)/%, where E is measured 
from the Fermi energy and 7 is a small imaginary part. 
Throughout this paper, we fix parameters as Ln = 70, 
W = 25, /j, = 2t, and Vj = 2t. Under these parame- 
ters, normal metals are in the diffusive transport regime, 
where the mean free path in normal metals is estimated 
about £ ~ 6 lattice constants and the Thouless energy 
Eth is calculated to be 1.6 x 10~ 3 t. Results discussed 
below are qualitatively insensitive to these parameters. 
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where cj, a (c r CT ) is the creation (annihilation) operator 
of an electron at r with spin a = ( f or J, ) and c means 
the transpose of c. The hopping integral t is considered 



FIG. 2: The maximum amplitudes of Josephson currents 
in the p^-wave symmetry J c (px) are compared with those in 
the s-wave symmetry J c (s) in (a), where Tb is the trans- 
mission probability of potential barriers in the normal states. 
In (b), J c {Px) and J c (s) are plotted as a function of Tb at 
T = 0.001T C . 

At first we show that the maximum amplitudes of 
Josephson currents in the p^-wave symmetry J c {p x ) be- 
come much larger than those in the s-wave J c (s). In 
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Fig. Ufa), ratios J c (p x )/Jc(s) are plotted as a function 
of temperatures for A = O.li. Here we choose sev- 
eral values of the barrier potentials Vb at j = 1 and 
Lpf. The resulting normal transmission probabilities of 
the barrier T B are 1.0, 0.075 and 0.013 for V B /t = 0, 
6 and 15, respectively. The ratios J c {px) / Jc(s) increase 
with decreasing T and amazingly become more than 100 
in low temperatures for small Tg. The amplitudes of 
Josephson currents in the p x -wave junctions are much 
larger than those in the s-wave junctions. In Fig. Hfb), 
J c Rn normalized by 7rA /e is plotted as a function of 
Tb at T = 0.001T C , where Rm is the normal resistance of 
junctions. The results show that J c Rn in the s-wave de- 
creases with decreasing Tb, whereas that in the p x -wa.ve 
increases. 
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FIG. 3: Current-phase relations for the p^-wave symme- 
tries are shown for several temperatures at Vb = 0, where 
Ao = O.Olt in (a) and Ao = O.OOOlt in (b). For comparison, 
results in the s-wave junctions at T — 0.001T C is shown with 
a solid line in (a), where the amplitude of Josephson current 
is multiplied by 5. 

We next focus on current-phase relations of the Joseph- 
son effect. In Fig. [31 Josephson currents are plotted as a 
function of ip for the p^-wave symmetries at Vb = 0. Pa- 
rameters are chosen as Ao = 0.01/j and 0.0001/j in (a) and 
(b), respectively. The current-phase relations are almost 
sinusoidal function in a high temperature at T = 0.5T C . 
At T — 0.001T C , however, the current-phase relations are 
close to J oc ip and J oc sin(<p/2) in (a) and (b), respec- 
tively. These are characteristic current-phase relations 
in ballistic Josephson junctions in the s-wave symmc- 
trjii2iiii4. We have confirmed that these current-phase 
relations remain even in the presence of potential barriers 
(i.e., V B ^0). _ 

The results imply large contributions of the multi- 
ple Andreev reflection in low temperatures. In general, 
Josephson currents can be decomposed into a series of 
J = Y^=i Jn sin(ncp), where J„ for n > 2 represent con- 
tributions of the multiple Andreev reflection. Roughly 
speaking, J n is proportional to {Tjv}™ with Tjv being 
the transmission probability of a quasiparticle from the 
left superconductor to the right superconductor through 



the normal segment (including two barriers and a normal 
metal) . Thus the multiple Andreev reflection is negligible 
(i.e., Ji > J 2 > J 3 . . .) for Tjy < 1. On the other hand 
in the case of Jjv = 1, the multiple Andreev reflection 
leads to the deviation of current-phase relations from the 
sinusoidal function. It is noted at Tn = 1 that we ob- 
tain J oc ip and J oc sin(cp/2) at the zero temperature for 
L*N >" Co and L N <C £0, respectivel y 12 ! 13 ! 14 . 

In Fig.[3fa), we also show the current-phase relations 
in the s-wave symmetry at T — 0.001T C with a solid line. 
The current-phase relation in the s-wave is described al- 
most by the sinusoidal function^ because impurity po- 
tentials in normal metals suppress Try and therefore the 
multiple Andreev reflection. In the p^-wave junctions, 
the coherence length £ are estimated about 50 lattice 
constants in (a) and 5000 in (b). Thus Lm > £0 and 
Ljv <C £0 are satisfied in (a) and (b), respectively. The 
current-phase relations such as J oc sin(co/2) in (b) and 
J oc cp in (a) are universal properties of the p^-wave junc- 
tions in low temperatures because they are independent 
of the strength of barrier potentials and the degree of dis- 
order in normal metals. The calculated results in Fig. 
indicate Tjv — 1 even in the presence of impurity poten- 
tials. The large amplitudes of the Josephson current in 
Fig. are also explained by Tn = 1 • 

The calculated results in Figs. |2| and show the spe- 
cific properties of Josephson currents in the p^-wave junc- 
tions. In what follows, we analyze quasiparticle states in 
normal metals to understand the origin of the anoma- 
lous Josephson effect. In Fig. 0] we show the local den- 
sity of states in normal metals for the s- and p^-wave 
symmetries, where Ao = 0.005i, 7 = 0.05Ao, and 
denotes the normal density of states. At ip — in the 
s-wave junctions in (a), the local density of states for 
E < E t h ~ 0.3Ao is suppressed because of the proximity 
effect. The suppression of the local density of states in- 
dicates the conversion of quasiparticles to Cooper pairs 
in normal metals. At ip — n in (b), the local density 
of states recovers its amplitude for E < E t h ■ The wave 
function of Cooper pairs from the left superconductor 
and that from the right one cancel each other around 
ip ~ 7r as schematically illustrated in a picture below the 
calculated results. 

The local density of states is drastically changed in the 
Parwave symmetry as shown in (c) and (d). Zero-energy 
peaks whose width is determined by 7 can be seen, which 
means formation of the midgap Andreev resonant state 
(MARS) in normal metals. Although the MARS origi- 
nally localizes at junction interfaces^, the MARS pene- 
trates into normal metals in the presence of the proxim- 
ity effect. Spatial profiles of the local density of states 
depend remarkably on the external phase difference as 
shown in (c) and (d). At cp = 0, the zero-energy peak 
disappears at the center of normal metals (j ~ 35) be- 
cause wave function of the MARS from the left super- 
conductors cancel out that from the right one as shown 
schematically in a lower pannel in (c). On the other hand 
in (d), wave functions of the MARS in the two supercon- 
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FIG. 4: Local density of states (LDOS) in normal metals 
(1 < j < Ln ~ 70) are shown for the s-wave and p^-wave 
symmetries. The left and right superconductors are attached 
at j = and j = 71, respectively. Noting that E t h is about 
0.3Ao. In schematic pictures, DNM and S denote a diffusive 
normal metal and a superconductor, respectively. The local 
density of states shown here are calculated in the absence of 
Josephson currents. We have confirmed that the results at 
ip — 0.99-7T qualitatively shows the same behavior as those at 

If = 7T. 



ductors have the same sign with each other. Thus the 
two MARS's can penetrate deeply into normal metals 
and the zero-energy peak can be seen everywhere. We 
note that the penetration of the MARS is possible only 
when the proximity effect is present in normal metals. In 
fact, we have confirmed that no zero-energy peak is found 
in normal metals in the d xy -wave symmetry (results are 
not shown) and that the ensemble average of Josephson 
currents vanishes because the proximity effect is absent 
in normal metals^. Fig. 0| indicates that the proximity ef- 
fect bridges the two MARS's in the two superconductors. 
Thus Tm — 1 holds because of the resonant transmission 
through the MARS in normal metals. The Josephson ef- 
fect specific to the p^-wave symmetry discussed in Figs.|2 
and|2]are a consequence of the diffusion of the MARS into 
normal metals. 

In summary, we found anomalous behaviors of Joseph- 
son currents in superconductor / normal metals / super- 
conductor junctions in the p x -wave symmetry. The max- 
imum amplitudes of Josephson currents J c in the p^-wave 
junctions become much larger than those in the s-wave 
junctions. It is known that large values of J c are de- 
sired in device applications because J c Rn limits opera- 
tion speeds of Josephson devices. Current-phase relations 
in low temperatures are close to those in ballistic junc- 
tions such as J oc sin(y/2) and J cx ip independent of the 
strength of potential barriers at interfaces and the degree 
of disorder in normal metals. The two the midgap An- 
dreev resonant states penetrate deeply into normal met- 
als, which causes the unusual Josephson effect in p^-wave 
superconducting junctions. The anomalous Josephson 
effect is a novel feature of phase-sensitive transport in 
spin-triplet superconducting junctions. 
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